Measurement of distance and orientation of two atoms in arbitrary geometry 
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Accurate measurement of relative distance and orientation of two nearby quantum particles is 
discussed. We are in particular interested in a realistic description requiring as little prior knowledge 
about the system as possible. Thus, unlike in previous studies, we consider the case of an arbitrary 
relative orientation of the two atoms. For this, we model the atom with complete Zeeman manifolds, 
and include parallel as well as orthogonal dipole-dipole couplings between all states of the two atoms. 
We find that it is possible to determine the distance of the two atoms independent of the orientation, 
as long as the particles are sufficiently close to each other. Next, we discuss how in addition the 
alignment of the atoms can be measured. For this, we focus on the two cases of atoms in a two- 
dimensional waveguide and of atoms on a surface. 

PACS numbers: 42.65.Sf, 42.50.Nm, 42.60.Da, 04.80.Nn 



I. INTRODUCTION 

Progress in many areas of science and its application 
is fueled by the ongoing progress to measure and struc- 
ture small objects. In many cases, light is used as a 
primary tool for reading or writing. But since light is 
subject to diffraction, a straightforward implementation 
is restricted to structures of order of the involved wave- 
length [H. Different methods have been invented to sur- 
pass this limit, such as near-field imaging 0, techniques 
based on the selective addressing of nearby particles ^SJ , 
resolution enhancement due to non-classical effects [J], 
multiphoton spectroscopy [F] , quantum lithography with 
classical fields [6] , or position-dependent dark states 01 ■ 

Among the most fundamental problems in this area 
is the measurement of the distance between two nearby 
quantum particles such as atoms. It has been recognized 
that a precise determination of the interparticle distance 
is possible down to distances far below the wavelength 
of the employed light based on their mutual interaction. 
For small distances, the atoms are coupled by the dipole- 
dipole interaction, which modifies the optical properties 
of the system Q . This was confirmed in a recent ex- 
periment [lo| . and it was found that the resonance fluo- 
rescence exhibits characteristic features which enable one 
to determine the relative distance over a large range of 
small distances The resonance fluorescence has the 
advantage that it can be observed in the far field, and 
distance determination via fluorescence is not affected 
by the usual resolution limitations since the distance in- 
formation is encoded in the frequency spectrum of the 
emitted light. Similar ideas have also been used for the 
localization of single particles [13 . 

The existing distance measurement techniques based 
on the dipole-dipole interaction, however, are restricted 
to two two-level atoms in specific geometries, such as 
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aligned along the propagation direction of the exciting 
laser field. In most practical cases, however, the relative 
orientation of the two nearby atoms is unknown, and for 
similar reasons, it is equally difficult to measure the rela- 
tive orientation as the distance. Thus the question arises, 
whether the ideas of [l0|, EH can be extended to the case 
of arbitrary orientation. It turns out that it is not mean- 
ingful to study the system of two two-level atoms in the 
case of arbitrary orientations. The reason for this is the 
appearance of dipole-dipole couplings between orthogo- 
nal transition dipole moments (DDOTDM) in more gen- 
eral geometries [H, [11, [H, [Tl, The electric field 
emitted by one of the particles has not only a compo- 
nent corresponding to the emitting transition dipole mo- 
ment, but also a component along the interparticle dis- 
tance vector. The projection of the latter field compo- 
nent on a transition dipole moment in the second atom 
can be non-zero even if it is orthogonal to the emitting 
dipole p^ . In a real atom with magnetic level struc- 
ture, these DDOTDM lead to the population of excited 
states even if they are not driven by the external laser 
field. Thus, the two-level approximation breaks down, 
and correct predictions can only be expected if the the- 
oretical modelling includes complete Zeeman manifolds 
including all occurring dipole dipole couplings [l^ . 

Motivated by this, here we study the determination of 
relative distance and orientation of two nearby atoms in 
arbitrary geometry. The atoms are driven by a single res- 
onant standing-wave laser field, and we make use of the 
far field resonance fluorescence intensity and spectrum 
as observables. Each atom is modelled as a four-level 
system with one ground state (total angular momentum 
zero) and three excited states (total angular momentum 
one), including all relevant dipole-dipole couplings oc- 
curring in arbitrary geometries. We start by analyzing 
the distance determination for the case of a known ori- 
entation, and present dressed-state interpretations of the 
obtained resonance fluorescence spectra in various cases 
of relevance. Next, we describe a method to determine 
the interparticle distance for arbitrary orientation, which 
works as long as the particles are sufficiently close to each 
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FIG. 1: (Color online) Scheme for the determination of rela- 
tive distance and orientation of two nearby atoms. The atoms 
H £ {1,2} are driven on one transition by a standing wave 
laser field with Rabi frequency fi(r,,) and scatter light, which 
is registered in the far field with a detector. The interatomic 
distance vector R is arbitrary, as shown in the left inset. The 
right inset shows the level structure of each atom. Each atom 
has a single ground state (zero angular momentum) and three 
excited states (angular momentum 1). 

other. Finally, we discuss methods to determine the rela- 
tive orientation of the two particles, focusing on the two 
cases of atoms confined in a planar waveguide and atoms 
on a surface. 

II. THEORETICAL CONSIDERATIONS 
A. Master equation 

Our system consists of two identical nearby atoms lo- 
cated at Ti {i G {1,2}) as shown in Fig. [1] Each atom 
is modelled as a four-level system with one ground and 
three excited states, modelling the complete magnetic 
substructure of a transition from an 5' state (angular 
momentum 0) to a P state (angular momentum 1). The 
Hamiltonian of the two atoms coupled to the surrounding 
vacuum field is given by 

H = Ha + Hp + Hyac , (1) 

where 

^^ = ^EE^'^1+^- ' (2a) 

Hp hukal^aks , (2b) 

ks 

H.ac = -d^'^ ■ E{ri) - d^'^ • E{r2) . (2c) 

In these equations. Ha is the Hamiltonian that describes 
the free evolution of the bi-atomic system. We set the en- 



ergy of the ground state to zero, and the energies of the 
excited atomic states \i) {i G {1, 2, 3}) are huJi. The rais- 
ing and lowering operators on the |4^) <-> transition 
of atom /i (/i S 1,2) are 

slt^ = M{U and ^i^^' = |4^)(v|. (3) 

The Hamiltonian of the vacuum field is described by Hp, 
with ttks and a^.^, as the field annihilation and creation 
operators. Hyac describes the interaction of the atom 
with the vacuum field in the dipole approximation with 
the vacuum field E{r) given by 

= ^ E + ^- ■ (4) 

ks 

k is the wave vector, eus the polarization, uj^ the fre- 
quency of a field mode, and V the quantization volume. 
We use the Wigner-Eckart theorem [18j to determine the 
electric dipole moment operator of atom fi which is given 
by 

d'^^=^d,5(^)+H.c., (5) 

i=l 

and the dipole moments di = {i\d\4:) are given by the 
matrix elements of the electric dipole moment operator 
d as 

di d2=Ve,, d3 = -Ve^-K (6) 

£(=•=) _ _|_ iey)/^/2, and V denotes the reduced dipole 
matrix element. Note that the dipole moments di are 
independent of the atomic index fi since the two atoms 
are identical. The vector that defines the relative position 
of the atoms in spherical coordinates is given by 

R = r2 — Ti = i?(sin 6 cos (j), sin 6 sin 0, cos 9)^ . (7) 

Using standard methods, the system Master equation 
evaluates to [H, [H, [H, ^ 

^^~^[HA,g]-^[Hn,g] + C^g. (8) 

The Hamiltonian H^^ arises from the coherent part of the 
dipole-dipole interaction and is given by 

3 

Hn = -hY^irtuS^l^ sll^ + H.c.} - h{n2i(s^2ls[^l 
1=1 

+ 5(^)5fJ) + H.c.} - Mf^3i(sfis(^J + S^sfJ) 
+ H.c.} - h{n,2iS^^lsi'l + S^^y^l) + H.c.} . (9) 
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The incoherent part of Eq. ^ is given by 
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+ Ts^iS^'-^si-^g + gS',^^S['2 - 2St!gS',^^ 

+ T,,isi^^si^_}g + gSi^^si^] - 2Si^_}gSi^^) + H.c.} . 

(10) 

In the above equations, coefficients fiy and Tij with i = 
j represent dipole-dipole couphngs between a transition 
dipole of one atom and the corresponding parallel dipole 
of the other atom, whereas terms with i j correspond 
to dipole-dipole couplings between orthogonal transition 
dipole moments in the two atoms. The coefficients fiy 
and Tij can be calculated as 
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= -[diRe{X)d*] 



(11a) 
(lib) 



Here, Re( X ) [Im( X )] represents the real [imaginary] 
part of the tensor X whose components X^i (fc, I G 
{1,2,3}) are given by 
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:^(i2) = i[4K- + 4-- 

47r£o '7 V V 
RkRi , 1 , 3t 3 

52~(~ + ~ 3)i' 



(12) 



6ki is the Kronecker delta symbol and 77 = koR, with the 
approximation wi « W2 ~ W3 « wq, where loq = ckg is the 
mean transition frequency. From Eq. (|llap . the different 
dipole-dipole coupling constants flij evaluate to 
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{[3r/2 - 1 + {rf - 3) cos 20] cos 77 



8773 

7](1 + 3cos26') sin 7/} 
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fl2i = —\/2 cot 6*1731 ( 

1^22 = rjii - (2 cot^ e - 1)^2316^"^ , 

^^32 — "^211 ^33 — ^11 • 
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FIG. 2: (Color online) (a) Magnitude of few orthogonal 
dipole-dipole coupling constants. The parameters are R = 
0.08A and = ■k/2. The red solid curve shows |f2i3|, while 
the blue dashed curve depicts |f22il ~ \^32\- (b) Magnitude 
of few parallel dipole-dipole coupling constants at S = tt. The 
red solid line shows 1^22!, and the dashed blue line depicts 

Ifilll = 1^^331. 



low from Eq. (jllbp as 
3 

47/ 



r3i = 7^ [(7/2 - 3) sin 7/ -I- 37/ cos 77] sin^ eie"^''^ , (14a) 



En ==3--^{[3772-l + (7/2_3)cos26i] sin 77 
87/ 

+ 77(1 + 3 cos 26*) cos 77}, (14b) 

E21 = -\/2cot6lE3ie"^, (14c) 

^22-^n-(2cot2^7-l)F3le2^^ (14d) 

r32 = -r2i, E33 = Fn. (14e) 

Some examples illustrating the dependence of the cou- 
pling constants flij on or i? are shown in Fig. [2] 

The total spontaneous decay rate of each individual 
atom is given by 27^, where 



7 := li 



1 2|d,pL^g 
47reo 3?ic-^ 



(15) 



and we have again used the approximation iOi « wo- 

We finally include an external driving laser field, which 
we assume to be polarized along the z-axis. Its electric 
field is given by 



The corresponding incoherent coupling constants Fy fol- 



(16) 
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where £z denotes the amphtude, ujl the frequency, and 
Ez is the polarization of the field, and c.c. denotes the 
complex conjugate. The wave vector k l with wave num- 
ber fci — 27r/A points along the positive x-axis. 
In a suitable interaction picture, we obtain 

^^=-'j[HL + HArQ]-^[Hn.g]+C^~Q, (17) 

with 

^^ = -^EE^»^l+^- • (18) 

/I — 1 2—1 

Here, the detunings /S.i — ujl — ivt {i e {1,2,3}). The 
interaction of the system with the external laser field in 
the electric dipole and the rotating wave approximation 
is described by 

HL = -^j2in{r,)si>:;^+R.c.]. (19) 

The position dependent Rabi frequencies are given by 

ri(r^) = r!sin(fcL -r^), (20) 

where ft = VEz/fi {fJ, £ {1,2}). Due to its polarization, 
the laser field couples only to the |2) ^ |4) transitions in 
the two atoms. 



B. Observables 

Our observables are the resonance fluorescence inten- 
sity and the resonance fluorescence spectrum of the light 
emitted by the two atoms. The spectrum of resonance 
fluorescence up to a geometrical factor is determined by 
real part of the Fourier transform of the two time corre- 
lation function of the electric field [l9l |. 

(21) 

In this equation, iE [E ] denotes the positive [nega- 
tive] frequency part of the electric field operator and the 
subscript st refers to the steady state. At a point r — rr 
in the far-field zone, the negative frequency part of the 
electric field operator evaluates to Q 

E^'\r,t) = E\rl{r,t) 

- i^T^ E'^' ^^(^^ S^+{iy"'\ (22) 

2 — 1 

where t — t— r/c is the retarded time and Si±(t) — 
exp{^LUJLt)Si±{t). The first term denoting the free field 
can be neglected if the point of observation lies outside 
the driving field. 



The resonance fluorescence intensity is given by the 
one-time normally ordered correlation function of the 
electric field operators, 

Ist = {E^'\r,t)-E^^\r,t)),t. (23) 

In the following, the point of observation is assumed to be 
along the y or —x direction for the resonance fiuorescence 
spectrum, and along the z direction for the resonance 
fluorescence intensity. Evaluating the cross products in 
Eq. (im, we find that our choice of observation direction 
enables us to separate linearly polarized light emitted on 
transitions |2) |4) from the circularly polarized light 
emitted on transitions \i) ^ |4) (i € {1, 3}) by means of a 
polarization analyzer. We designate linearly [circularly] 
polarized spectra as tt [a] ones, respectively. 

III. RESULTS 

The resonance fluorescence spectrum Eq. (I?T|) emitted 
by the two atoms in general is rather complicated, but 
it simplifies considerably in certain parameter cases, as 
it was found already in the case of two nearby two-level 
systems [ll| . In the following, we will in particular re- 
fer to the case of either small or large interatomic spac- 
ing, on a length scale given by the involved transition 
wavelength. For small distance, the coherent part of the 
dipole-dipole interaction dominates the system dynam- 
ics, with corrections due to the much weaker laser field 
Rabi frequencies. In the opposite case of larger sepa- 
ration, the Rabi frequencies dominate, with corrections 
from the dipole-dipole interaction. Spectra for situations 
in which the dipole-dipole interaction and the Rabi fre- 
quencies are comparable usually can not be interpreted 
in a straightforward way. In these cases, the driving field 
intensity can be increased or decreased in order to evolve 
in one of the two simpler cases. In the following, we will 
make use of this general observation, and present our re- 
sults in two steps. First, we will describe methods to 
determine the interatomic separation in various cases of 
interest. Second, we will discuss the determination of the 
relative orientation of the two atoms. 

For the numerical analysis, we assume that ri — 
(0.05 A, 0, 0). Our measurement techniques, however, also 
apply to other values of ri. A special case arises if one of 
the atoms is at a node of the standing wave field. Such 
situations can be circumvented by shifting the phase of 
the standing wave slightly. We also assume the resonance 
condition, i.e., Ai = 0, see Eq. (fT5)) . 

A. Determination of the interatomic separation 

1. Known orientation of the atoms 

As a first step, we will present results for the situation 
of a fixed, known orientation. In particular, for simplic- 
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FIG. 3: (Color online) (a) Incoherent spectrum of resonance 
fluorescence for the case of larger atomic distances. The pa- 
rameters are R = 0.3A, (/) = 0, = 7r/2 and H = IOO7. The 
side peaks are located at the Rabi frequencies, f2(ri) « 30.907 
and r2(r2) ~ 8O.9O7. (b) The dressed state representation of 
the system in (a) showing the ac-Stark splitting of the two 
atoms. The splitting in each atom corresponds to the respec- 
tive position dependent Rabi frequency. 



ity, we analyze the case 9 = 7r/2,0 = 0. Then, the or- 
thogonal dipole-dipole couplings il2i and {i G {1,3}) 
vanish. Thus, the population is trapped only in the lev- 
els 1 2^) and |4^), and the system essentially reduces to 
that of two two-level atoms fu^ . This simple case, how- 
ever, suffices for the main aim of this section, which is to 
extend previous results by interpreting each case in de- 
tail in terms of the corresponding dressed state picture. 
This insight will enable us to explain our later main re- 
sults for general geometries. In principle, the results in 
this section also generalize to more complicated known 
orientations. While then the dressed-state analysis is in 
complete analogy to our discussions here, the analytical 
expressions are considerably more complicated. In any 
case, if the orientation is known, a numerical fit of the 
measured spectrum leads to the desired distance infor- 
mation. In the following Section fill A 21 we will extend 
our analysis to arbitrary orientations. 

For larger interatomic separation (about A/ 10 < R < 
A/2), the dipole dipole interaction is almost negligible 
and the two atoms can be considered independent. A 
typical spectrum is shown in Fig.[3Ua), and consists of five 
peaks. The middle peak occurs at w^, and the two side 
peak doublets are symmetrically located at uj^ ^^irf^)- 
An interpretation of the peak structure in terms of the 
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FIG. 4; (Color online) (a) Incoherent spectrum of resonance 
fluorescence for intermediate distances. The parameters are 
R = 0.08A, 6^ = 7r/2, = and f2 = 2OO7, such that the laser 
Rabi frequencies are larger than the relevant dipole-dipole 
coupling constants. The centre of the side band structures are 
located at the respective position dependent Rabi frequencies, 
±Q.{rx) « ±61.87 and ±r2(r2) « ±145. 79A. (b) Dressed-state 
representation, (i) shows the system without couplings, (ii) 
includes the dominant laser field couplings, (iii) in addition 
includes corrections due to the dipole-dipole interaction. 



system dressed states is shown in Fig.[3Jb). It turns out 
that the sidebands are shifted by the position-dependent 
Rabi frequencies experienced by the two atoms, respec- 
tively. Thus, the peak positions give a direct measure of 
the position of each atom and hence the interatomic sep- 
aration can be calculated using the expressions for the 
position-dependent Rabi frequencies Eq. ([^U)) pl| . The 
atomic levels in atom /.t split by the amount ri(r^), just 
as in the MoUow spectrum [20|]. Note that the two split 
doublets may coincide, for example, for 6 G {G,7r}, since 
then Vl{ri) = f^(r2). 

For intermediate interatomic distances (approx. 
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FIG. 5: (Color online) Incoherent spectrum of resonance flu- 
orescence in the small distance case, with laser field Rabi 
frequencies negligible compared to the dipole-dipole coupling 
constants. The parameters are R = 0.04A, 9 = n/2, (j> = 
and n = 2O7. The side peaks occur at tJL ± ^22- (b) Dressed 
state representation. The energies of the symmetric |s) and 
the antisymmetric \a) states are shifted by the dipole dipole 
coupling parameter 0,22- 



A/30 < R< A/10), the dipole dipole interaction becomes 
relevant. If the parameters are such that the driving 
field Rabi frequency is comparable to the dipole-dipole 
interaction, then the spectrum obtained is difficult to in- 
terpret. This, for example, occurs if the distance is de- 
creased to R = 0.08A in Fig. [3)^ a). We thus increase the 
Rabi frequency of the driving laser field and obtain the 
spectrum shown in Fig.lDJa). It can be observed that the 
sidebands are split by the dipole-dipole interaction. The 
corresponding dressed-state picture is shown in Fig.^Jb). 
This picture is obtained by evaluating the eigenstatcs and 
eigenvalues of the interaction picture Hamiltonian in the 
limit n{ri) ^ ^22- We consider the case il{ri) ^ fl{r2) 
and assume ft{r2) > 17(ri). In part (i) of Fig. Ufb), there 
is no interaction between the two atoms, such that four 
atom-field states are degenerate. In (ii), we include the 
interaction with the strong driving field, which splits the 
level scheme into four dressed states. The correspond- 
ing dressed states are listed in Tab. HI Finally, in part 
(iii), the weak dipole-dipole interaction is added, which 
further shifts the eigenenergies of all dressed states by 
approximately ±1722/2. 

There are sixteen possible transitions between the 
eigenkets having n + 1 and n photons. Transitions 
\p,n+l) \p,n), |;,n+l) ^ and \q,n + l) \q,n) 



\p,n + l) 
\m,n + 1) 
\q,n + l) 
\l,n + l) 



(1, 1, 1, l)/2 ^(^(ri) + Q{r2) + ^22)/2 
(1, -1, -1, l)/2 -^(^(ri) + Q.{y2) - ^22)11 
(-1, -1, 1, l)/2 ri(n(ri) - !^(r2) - !^22)/2 
(-1, 1, -1, l)/2 -^(^(ri) - n(r2) + 0.22)11 



TABLE I: Eigenvectors and eigenvalues of the interaction 
Hamiltonian of two two- level atoms in the limit f2(ri) ^ 0.22- 



and |m, n + 1) ^ |m, n) correspond to transition fre- 
quency wl. The eight transitions having frequencies 
± (f7(ri) ± ^22), i G (Ij^) form the two prominent 
side band doublets on each side (eight side peaks in to- 
tal) visible in Fig. [^a). These peaks are crucial for the 
distance determination. Setting parameters such that 
fi(ri) = 17(r2) would reduce the number of transitions 
± (ll(rj) ± ^22) to four. The remaining four transi- 
tions involving both Rabi frequencies WL±(ri(r2)±ll(ri)) 
are hardly visible. 

In this case, distance measurement is possible in two 
ways. First, the center of the sideband peak doublets 
correspond to 17 (r^), such that again a direct position 
determination of the two atoms is possible via Eq. (|20p . 
Second, the doublets are split by the dipole-dipole cou- 
pling strength ^22- Thus, from Eq. (jl3d[) . again the dis- 
tance can be obtained. Best results are obtained by a 
combination of the two methods. 

At small interatomic separation (approx. i? < A/30), 
the dipole-dipole interaction typically dominates the sys- 
tem dynamics. In Fig.[S]Ja), two peaks can be observed in 
the spectrum of resonance fluorescence which occur at the 
frequencies cj^ ± ^22- These peaks reflect the well-known 
symmetric and anti-symmetric states (\eg) ± \ge))l\f2 
formed in a dipole-dipole interacting system of two two- 
level atoms as shown in Fig. [D^b) . The two peaks corre- 
spond to the transitions |e) — > |s) and |s) — > while 
the anti-symmetric state decouples from the dynamics in 
the limit i? ^ 0. Thus, from the positions of the peaks, 
one can immediately obtain the dipole-dipole interaction 
energy and thus the interatomic distance via Eq. (flSdll . 

We finally analyze the case in which a laser field per- 
turbs the spectrum generated due to a strong dipole- 
dipole interaction. In Fig. [DJa), the driving laser field 
with Rabi frequency at the anti-nodes 17 = 757 leads to 
17(ri) = 23.I87 and 17(r2) = 40.197, which is smaller 
than the magnitude of relevant dipole-dipole coupling 
1722 = 91.647. Instead of two side peaks, one now ob- 
tains split peaks at the sides as well as in the middle. 
The dressed state picture of this situation is shown in 
Fig. mb), assuming 17(r^) ^ 1722- Starting from the 
non-interacting system in (i), in part (ii), the dominat- 
ing dipole-dipole interaction is included. It combines 
the states leg, 71 + 1) and |(?e,n + 1) to form |s,n + 1) 
and |a,rt + 1) as symmetric and anti symmetric com- 
binations, respectively, and shifts their energies by the 
amount of the dipole-dipole interaction energy. In part 
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FIG. 6: (Color online) (a) Incoherent spectrum of resonance 
fluorescence. The parameters are as in Fig. [5] except for 
Q, = 757, such that the dynamics is dominated by the dipole- 
dipole interaction, but notably perturbed by the driving field, 
(b) Dressed state representation, (i) uncoupled states, (ii) in- 
cludes the dipole-dipole splitting, (iii) indicates the additional 
coupling to the laser field, and (iv) shows the full dressed 
states induced by the dipole-dipole coupling perturbed by the 
laser field Rabi frequencies, £7i means n(ri). 



(iii), the laser field is included into the dynamics. Af- 
ter another basis transformation to also dress the system 
with the laser field, this results in the further shifting of 
the eigenstates by approximately [VI{y i) ±^(t / {^^22), 
as shown in (iv). The corresponding eigenenergies and 
eigenstates are listed in Tab. |lTl 

Sixteen transitions between the eigenkets having n + 1 
and n quanta take place. Transitions |0,n -I- 1) — > 
|0,n), |+,n + l) \+,n), \-,n + I) -> |-,n) and 
|a,n -t- 1) ^ occur at transition frequency uol- 

Transitions |0, n+1) 1 + ,?^) and | + ,n+l) — > |0,?t-) 
have respective frequencies ± ^22 ± i^(ri)ri(r2)/ri22- 
The transitions |0,n -I- 1) ^ | — j") and | — ,n -I- 1) ^ 
|0,n) involve frequency differences equal to ujl ± ^22 ± 
(ri(ri) -I- ri(r2))^/(2Sl22), respectively. The separation 
in frequency for the transitions |0,n-f 1) — > \a,n) and 
|a, n-l- 1) |0, n) approximates to ±2Jl22 ± (il^(ri) + 
f^^(r2))/(2ri22), respectively. Finally the transitions 



State Composition Energy 

la, n + l) (0, -1, 1, 0)772 -h{Q,22 + (n(ri) - n{v2)f / {^22)) 

| + , n + l) (-1, 0, 0, 1)772 ft(n(ri) - n{v2)f / {^^22) 

h,n + l) (1,0,0,1)772 -/i(n(ri) + n(r2))V(4n22) 

|0,n + l) (0,1,1,0)772 fi(n22 + (»(ri)+»(r2))V(4n22)) 



TABLE II: Eigenvectors and eigenvalues of two identical two- 
level atoms in the limit f2(r^) <^ 0.22- 



| + ,n+l) \ — ,n) and | — ,n-f 1) correspond to 

the frequency difference ± (f2^(ri) -|- f2^(r2))7(2f222), 
respectively. These transition frequencies are the posi- 
tions of the peaks in the resonance fluorescence spectrum 
of Fig.[nija). The corresponding frequencies for the tran- 
sitions | + ,n + l) \a,,n) and \a,n + 1) are 
bJL ± ^^22 ± (f^Cri) - ^{T2)f / (2^22) and for the tran- 
sitions | — ,n-|-l) \<i,n) and |a,n-|-l) \ — are 
ojl ± f^22 T f^(ri)f^(r2)7^22, respectively. The last four 
transitions do not show up in the spectrum. 



2. Unknown orientation of the atoms 

The methods presented so far allow for a determi- 
nation of the interatomic distance if the orientation of 
the two atoms is known and fixed. Often, however, the 
orientation is unknown. Therefore, in this section, we 
turn to our main results, and present a method to ob- 
tain the interatomic distance for arbitrary orientations. 
This method can be applied if the interatomic distance 
is sufficiently small, such that the dipole-dipole interac- 
tion dominates the system dynamics. Fortunately, this 
usually is exactly the parameter range in which a dis- 
tance determination is desired. In order to explain the 
method, we first imagine the two atoms without any driv- 
ing fields. Then, it turns out that the eigenenergies of the 
dressed states are independent of the orientation of the 
two atoms [14|] . An interpretation of this fact is that with- 
out external field, there is no preferred direction in space, 
such that the energies cannot depend on the orientation. 
Since the spontaneously emitted light is emitted at fre- 
quencies corresponding to the dressed state energies, it 
follows that the positions of the peaks in the fluorescence 
spectrum of the atoms are unaffected by the orientation 
of the atoms. This property is approximately preserved if 
the atoms are driven by a weak driving field, which has 
a Rabi frequency much smaller than the dipole-dipole 
couplings. We therefore find that at small distance and 
weak driving, the resonance fluorescence spectrum has 
peak positions independent of the alignment of the two 
atoms. Only the relative widths and heights of the spec- 
tral features change with the orientation. It is important 
to note that these properties of the two-atom system are 
only described correctly if all dipole-dipole couplings are 
included in the modelling [14]. This is the reason why we 
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FIG. 7: (Color online) Incoherent spectrum of resonance flu- 
orescence for small distance and weak driving fields. The 
parameters are K = 0.04A, <^ = 7r/15, Q — vr/S and = 2O7. 



included complete Zeeman manifolds in our analysis. 

It remains to deduce the interatomic distance from the 
peak positions in the resonance fluorescence spectrum. 
For this, we again analyze the eigenvalues of the inter- 
action Hamiltonian, which determine the peak positions. 
Since the peak positions and thus these eigenvalues are 
independent of the orientation, it suffices to evaluate an- 
alytic expressions for the peak positions in a simple con- 
figuration. Investigating the eigenvalues for 9 = 0, and 
assuming ri(T'^) <^ flu, we find that the eigenenergies 
are given by 0, ±ilii{6 = 0), and zLil22{S ~ 0). An ex- 
ample for this is shown in Fig. [71 The four side peaks are 
located at iiln and ±^22, and the interatomic distance 
can be gained from both coupling constants via Eq. (|13bp 
and (|13dp . As can be seen from Fig. [21 for 6* e {0, tt}, 
the coupling constant |r222| is larger than jrinj for small 
interatomic distances R. Thus, the inner [outer] peaks in 
Fig. [7] correspond to |riii| [|ri22|]- Due to the dependence 
of the amplitude of the spectral peaks on the orientation, 
the peaks at ztrin visible in Fig. [7] may be suppressed. 
For example, in Fig. [5ja), we found only a single pair of 
sidebands corresponding to ±D,22- This can be under- 
stood by observing that these peaks at ±il22 correspond 
to states 1 2) populated by the driving laser field, while 
states |1) and |3) are only populated in certain geome- 
tries. 

We end this section by noting that in the case of large 
interatomic separation, the dipole-dipolc interaction van- 
ishes, such that the interatomic coupling also becomes 
independent of the orientation. The driving field Rabi 
frequencies experienced by the two atoms, however, may 
not be the same, as they depend on the scalar product 
Icl • r^. Thus, by applying a driving field with larger 
Rabi frequency, it is possible to measure the position of 
the atoms projected on the propagation axis of the driv- 
ing field. For arbitrary orientations, however, a single 
measurement of this type does not allow to deduce the 
interatomic distance since the position transverse to the 
wave vector remains unknown. 

We thus conclude that by applying weak driving fields, 




detector 

FIG. 8: (Color online) Setup with the two nearby atoms con- 
fined inside a two-dimensional waveguide. The atoms are lo- 
cated in the x-y plane. 



the interatomic distance can be measured from a pair of 
nearby atoms independent of their mutual orientation, as 
long as the dipole-dipole interaction is strong enough to 
dominate the system dynamics. 



B. Determination of the orientation 

So far, we have discussed techniques for the measure- 
ment of the interatomic distance, and have demonstrated 
how the distance can be measured independent of the 
orientation of the interparticle distance vector. In this 
Section, we augment our analysis by discussing the deter- 
mination of the relative orientation of the two atoms. We 
discuss two different cases, corresponding to two different 
methods to determine the orientation. First, we discuss 
the case of unknown (p, assuming 9 = 7r/2. This case 
corresponds to an effective two-dimensional geometry of 
the system which can be realized, e.g., hy embedding the 
atoms in a planar matter waveguide [Hf, [H, [2^ . In this 
case, the orientation is deduced from the (/)-dependent 
peak positions in the fluorescence spectrum induced by 
the driving laser field. Second, we study the case of un- 
known 6 and (j) = 7r/2. This corresponds to atoms on a 
surface, driven by a laser field propagating perpendicular 
to the surface. In this case, we will determine 9 via the 
resonance fluorescence intensity emitted in a particular 
direction. 



1. Unknown <j): Planar waveguide 

In this section, we assume that the two atoms are con- 
fined in the x-y plane (6* = 7r/2) as shown in Fig. [51 as it is 
the case, for example, in a planar waveguide. Since such 
a setup may also constrain the observation direction, we 
assume detection in a direction anti-parallel to the inci- 
dent driving field, i.e., along the —x direction. This way, 
the spectrum can be measured without background from 
the incident laser field. In this geometry, the coupling 
constants fl2i and O32 are zero for all values of cj). Since 
only the second transition |2) ^ |4) is driven, the pop- 
ulations of the levels |1) and |3) are zero for the whole 
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FIG. 9: (Color online) Incoherent spectrum of resonance fluo- 
rescence for two atoms in a geometry as shown in Fig. [S] The 
parameters are R — 0.07A, 9 = n/2, (f> — O.lvr, and fl = 35O7. 



range of (f>. The parallel coupling constants flu = ^l^s 
and ri22 are independent of (f>. Nevertheless, the spectra 
depend strongly on <j) because of the ^-dependence of the 
Rabi frequency 51 (r2). The obtained spectra are identical 
if (j) is replaced by 27r — 0. 

From the results of Sec. lIII A2l it is clear that the peak 
positions in the resonance fluorescence spectrum cannot 
be used to determine the orientation as long as the dipole- 
dipole interaction dominates the dynamics. Therefore, 
we apply stronger driving fields, such that the external 
driving dominates the dynamics. 

In this case, for most values of 4>, a typical spectrum ob- 
tained is shown in Fig. [SI Using the results of Sec. lIII ATI 
we conclude that the two doublets on each side corre- 
sponding to MoUow sidebands at the two Rabi frequen- 
cies experienced by the atoms, split by the dipolc-dipole 
interaction. The doublets can thus be used to approx- 
imately read off the two position-dependent Rabi fre- 
quencies. Assuming that the distance is known from a 
measurement with a weaker driving field as described in 
Sec. nil A 21 the components of the position vectors of 
the individual atoms as well as the relative alignment of 
the atoms along the laser can be found out by using the 
position dependent Rabi frequencies ri(r^) as follows; 



1 



-1 f ^jr,] 



(24) 

In Fig. [9l the peaks in the spectrum of resonance flu- 
orescence occur at approximately 0, ±91.687, ±123.957, 
±218.597, and ±250.867. Considering the mean values 
of the peak separation in the inner side band doublets 
as f2(ri) and that in the exterior side band doublets as 
r2(r2), from Eq. we obtain 0.09l7r as the value of (p, 
which deviates from the true value by about 9%. This 
deviation can be attributed to the imperfect determina- 
tion of the Rabi frequencies as the mean value of the two 
peaks in the doublets. 

Increasing from the value that has been used in 
Fig. [21 the position-dependent Rabi frequencies change. 
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FIG. 10: (Color online) Fluorescence intensity emitted on the 
cr-transitions observed in z direction from two atoms on a sur- 
face. A9 is the relative angle between laser polarization and 
interatomic distance vector. The green dashed lines indicate 
the four possible values of AO corresponding to a possible 
measured intensity indicated by the solid red line. The pa- 
rameters are R — 0.07A, (j) = -n /2 and Q. = 2OO7. 



until the sideband doublets corresponding to the two 
Rabi frequencies start to overlap, as the two position- 
dependent Rabi frequencies approach each other. In this 
case, it is difficult to estimate cj) directly from the spec- 
trum, since a clear identification of the dilferent peaks 
is not obvious. One strategy is to increase the driving 
field intensity. Since the peak separation |ri(ri) —Q,{y2)\ 
is proportional to 17, this increase eventually leads to a 
splitting larger than the line widths of the involved peaks, 
such that an identification becomes possible. In any case, 
it can be concluded from overlapping peaks that (j) is close 
to 7r/2 or 37r/2, since then ri(ri) « ^{v2)- 

In summary, the relative orientation of the two atoms 
with respect to the laser can be determined using stronger 
laser fields. This works well if the position dependent 
Rabi frequencies are different from each other, since then 
the spectral lines are well separated. This is the case for 
not close to 7r/2 or 37r/2. Accordingly, if the correspond- 
ing spectral peaks overlap, it can be concluded that (j) is 
close to 7r/2 or 37r/2. 



2. Unknown 9: Atoms on a surface 

In this section, we consider the case of two atoms on 
a surface, driven by a laser field propagating perpendic- 
ular to the surface. Thus, is fixed to 7r/2, while is 
unknown. In this case, the two Rabi frequencies expe- 
rienced by the atoms are equal and independent of R 
and 0. Therefore, in contrast to the previous Sec. IIII B 11 
here, we determine the value of 9 with the help of the 
resonance fluorescence intensity. In particular, we con- 
sider the a-intensity emitted by the dipoles di and ^3, 
measured by a detector placed in z direction since there 
is no cr-spectrum in y direction. 

It turns out that the configuration is symmetric in the 
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sense that a rotation of the laser polarization and the 
detectors around the x direction is equivalent to a cor- 
responding rotation of the interatomic distance vector. 
Therefore, the measured resonance fluorescence inten- 
sity depends only on the relative angle A9 between the 
laser polarization direction and the orientation of the two 
atoms 6 on the surface. In Fig. (TUl we show this reso- 
nance fluorescence intensity versus the relative angle A6. 
A plot like this can either be recorded by rotating the 
sample in the y-z plane, or by rotating the polarization 
vector of the laser field around its propagation axis. 

From Fig. 1101 we find that the intensity is symmetric 
around A9 — 7r/2, and it is easy to see that one value of a 
intensity corresponds to at most four values of A9. These 
four values of can be roughly divided into the four 
ranges 7r/4, n/i 7r/2, 7r/2 Sn/i and Sn/i 
TT, respectively. This can be understood by noting that 
the orthogonal coupling constants responsible for the a 
intensity are the same for orientations 9, tt/2 — 9, tt/2 + 9 
and 11 — 9. At A9 e {0,7r/2,7r}, the a intensity is zero 
since the orthogonal dipole dipole coupling constants 
and O21 vanish at these points, see Fig.[21[a). Thus there 
is no population in states |1) and |3), and the intensity of 
light emitted from these states is zero. The points of zero 
a intensity A9 e {0,7r/2,7r} correspond to situations in 
which the polarization vector is parallel, perpendicular, 
or anti-parallel to J?, respectively. Since these values 
of A9 with vanishing intensity can easily be identified, 
they allow to determine 9 from the amount of sample or 
driving field polarization rotation required to reach these 
values. In particular, the symmetry point A9 ~ 7r/2 is 
well-suited for such a measurement. 



IV. SUMMARY AND DISCUSSION 

We have discussed methods to measure the relative dis- 
tance and orientation of two nearby atoms in arbitrary 
geometry. Our methods are based on the driving of the 
two atoms with a standing wave field, and on detection 
of the resonance fluorescence intensity and spectrum in 
the far field. The distance and orientation information is 
encoded in the scattered light via the position-dependent 
Rabi frequencies and via the distance- and orientation- 
dependent dipole-dipole couplings. Since unlike in pre- 
vious studies, we consider the case of arbitrary orienta- 
tion, the atoms must be described using complete Zee- 
man manifolds in order to correctly model all relevant 
dipole-dipole couplings between parallel dipole moments 
as well as between orthogonal ones. 

As preliminary work, we have analyzed the fluores- 
cence spectra in particular known geometries, in order 
to identify dressed-state interpretations in the various 
limiting cases of interest. These in particular are the 
case of dominating laser-induced dynamics perturbed by 
the dipole-dipole interaction, and the case of dominat- 



ing dipole-dipole coupling modified by the presence of a 
weaker laser field. Next, we have shown that the case 
of dominating dipole-dipole interaction enables one to 
measure the distance between two nearby particles inde- 
pendent of the relative orientation. The reason for this 
is that the eigenvalues of the total Hamiltonian describ- 
ing the dynamics, and thus the position of the system 
dressed states, are independent of the orientation if the 
two atoms are undriven. We found that a weak driving 
field allows to probe these dressed states without per- 
turbing the independence on the orientation. Finally, we 
discussed the measurement of the relative orientation. 
We presented two methods. The first is based on the 
position-dependent Rabi frequencies, which under cer- 
tain conditions reveal the orientation of the two particles. 
The second method is based on the measurement of the 
resonance fluorescence intensity in a particular direction. 
This intensity is a measure for the population in the ex- 
cited states not driven by the laser field, and therefore 
a signature for the magnitude of the dipole-dipole cou- 
pling between orthogonal dipole moments. We applied 
the two methods to the two cases of atoms confined in a 
two-dimensional waveguide, and to atoms on a surface, 
in which either the polar or the azimuthal angle of the 
interatomic distance vector is known. 

In principle, these methods to determine the orienta- 
tion can also be applied for the determination of both po- 
lar and azimuthal angle. The most promising ansatz is to 
make use of the resonance fluorescence intensity in a par- 
ticular direction as discussed in Sec. IIII ITS] together with 
a rotation of the sample or of the driving laser polariza- 
tion in order to fix one of the two angles at a value which 
renders the spectrum simpler (e.g., 0, 7r/2 or tt). Then, 
the methods described in Sees. IIII B II and l lll B 21 can be 
used to determine the other angle. The most straight- 
forward implementation, however, strongly depends on 
the experimental possibilities to modify the setup. For 
example, in many cases, a rotation of the sample will be 
difficult. 

For many applications, the generalization to more than 
two particles is desirable. It remains to be seen whether 
methods based on the dipole-dipole interaction can also 
be applied in such cases. One approach could be to com- 
bine methods presented here together with a selective ad- 
dressing of individual atoms at least in one or two dimen- 
sions, for example, by position-dependent state transfer. 
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